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ABSTRACT 

Optical photometry is presented for the quadruple gravitational lens 
PG1115+080. A preliminary reduction of data taken from November 1995 to 
June 1996 gives component "C" leading component "B" by 23.7 ±3.4 days and 
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components "Al" and "A2" by 9.4 days. A range of models has been fit to 
the image positions, none of which gives an adequate fit. The best fitting and 
most physically plausible of these, taking the lensing galaxy and the associated 
group of galaxies to be singular isothermal spheres, gives a Hubble constant of 
42 km/s/Mpc for Q = 1, with an observational uncertainty of 14%, as computed 
from the B — C time delay measurement. Taking the lensing galaxy to have 
an approximately E5 isothermal mass distribution yields H = 64 km/sec/Mpc 
while taking the galaxy to be a point mass gives H = 84 km/sec/Mpc. The 
former gives a particularly bad fit to the position of the lensing galaxy, while 
the latter is inconsistent with measurements of nearby galaxy rotation curves. 
Constraints on these and other possible models are expected to improve with 
planned HST observations. 

Subject headings: cosmology: distance scale, gravitational lensing — quasars, 
photometry 

1. Introduction 

Among the most promising applications of gravitational lenses is their use in 
determining cosmological distances (Refsdal 1964; Blandford and Narayan 1992). One 
measures fluxes for the components of a multiply imaged variable object as a function of 
time and determines the differential time delays. Though straightforward in principle, both 
the measurement of time delays and their interpretation has proven difficult in practice. 
Only recently have measurements for 0957+561 (Haarsma et al. 1996a,b and Kundic et al. 
1996 confirming the result of Schild and Cholfin 1986) and their interpretation (e.g. Grogin 
and Narayan 1996) begun to converge. 

The quadruply imaged z = 1.722 quasar PG1115+080 (Weymann et al. 1980; 
see Figure 1) has many properties which recommend it as a candidate for time delay 
measurements (Schechter 1996). But while PG1115+080 is known to vary (Vanderriest et al. 
1986), and though the flux ratios of the components appear to have varied (Schechter 1996), 
as yet no time delay has been reported. In the following sections we report measurement 
of multiple time delays from a preliminary reduction of data taken from November 1995 to 
June 1996. 
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2. Observations and Reduction 

Direct V filter CCD exposures were taken with the Hiltner 2.4-m telescope, the WIYN 
3.5-m telescope, the NOT 2.5-m telescope and the Dupont 2.5-m telescope, always in 
multiples of at least three exposures, typically 6 minutes each. Frames were bias subtracted 
and flatfielded and charged particle events were identified and flagged. Six bright stars were 
used to obtain the image scale and rotation. 

Photometry was carried out using methods nearly identical to those used by Schechter 
and Moore (1993) for the lens MG0414+0534. In the present case positions for the lensing 
galaxy and the four quasar images relative to each other were established from V exposures 
taken in 1994 with 0'/75 seeing, as were the shape parameters and the galaxy flux. The 
galaxy shape was convolved with the seeing appropriate to each image. Fixing the relative 
positions and the galaxy shape and flux leaves six parameters to be determined from each 
frame - the four quasar fluxes and the overall position of the system. 

Quasar fluxes were obtained relative to the nearby star designated "C" by Vanderriest 
et al. (1986), which was used as an empirical point spread function (PSF) template. Fluxes 
for other nearby stars relative to star "C" were also measured. No attempt was made to 
correct for color dependent airmass terms. 

Fluxes were summed for the two brightest components, Al and A2, which are separated 
by only Of! 48. Little is lost because the time delay between Al and A2 is expected to be 
hours. Fluxes for each multiplet of exposures were converted to relative magnitudes and 
averaged together, and an rms scatter computed. The median scatters for components 
Al + A2, C and B were, respectively, 4, 7 and 13 millimagnitudes. These were roughly 50% 
bigger than the formal errors in the fits. 

The average magnitudes for the quasar components and for star "B" are plotted^] in 
Figure 2. They include 30 nights of data from MDM Observatory and 2 from the NOT. 
Multiple entries on the same night indicate exposures taken with slightly different "V" 
filters. 

Fifteen nights of data taken with the WIYN 3.5-m telescope were also reduced but are 
not included. The seeing at the WIYN and the intranight consistency of these data were 
superb. However the night-to-night consistency was poor, and it has not proven possible 
to bring the MDM data and the WIYN data into accord. We expect that by allowing for 
variation in the PSF across the chip we will ultimately be able to combine the two data 



17 These data can be found in tabular form at 
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sets. 



3. Photometric Results 

In the light curves for components Al + A2, C and B shown in Figure 2, the errors in 
the means, as computed from the intranight scatter, are not much larger than the symbols. 
The night-to-night agreement is considerably worse; likewise the agreement of the shifted 
lightcurves is much worse than expected from the intranight scatter. On at least some 
nights the fluctuations, and those of star "B", appear to be correlated (e.g. JD 24450088). 
It seems likely that these are due to errors in the adopted PSF. Therefore in plotting the 
quasar components in Figures 2 and 3, and in estimating time delays, we have subtracted 
one half the magnitude of star "B" (taken with respect to star "C") from our magnitudes 
for the quasar components (also taken with respect to star "C"). This gives equal weight to 
stars "B" and "C" as photometric references. 

The amplitude for components Al + A2 appears somewhat smaller than that for C 
or B, but this is not unexpected. Witt, Mao and Schechter (1995) find that unresolved 
components suffer rms fluctuations of 0.6 mag due to microlensing in quadruply imaged 
systems. If the time variable core is unresolved by the microlenses, while most of the flux 
arises from a larger, resolved region, one expects different fractional flux changes in the 
different components. 

4. Time Delay 

We have computed time delays and flux ratios for the (Al + A2) — C and B — C pairs, 
using the method of Press, Rybicki, and Hewitt (1992), but generalized so that the four 
parameters are fit to the three light curves simultaneously. Since the errors derived from the 
intranight scatters are clearly too small, we have added to these, in quadrature, additional 
errors of 4, 7, and 10 millimag for components Al + A2, C, and B respectively. The present 
analysis incorporates the implicit assumptions a) that the errors are uncorrelated and b) 
that the fractional flux variations are the same for each component. 

Under these assumptions, we find component C leading component Al + A2 by 9.4 
days, and component C leading component B by 23.7 days. The best-fit C/(A1 + A2) 
flux ratio is 0.153 and the best-fit B/C flux ratio is 0.634. It follows that component 
Al + A2 leads B by 14.3 days and that the B/(A1 + A2) flux ratio is 0.097. Figure 3 shows 
lightcurves for components B and C, with the latter delayed. 
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A plot of the x 2 statistic for the (Al + A2) — C and B — C delays shows a clean 
elliptical minimum, indicating an error for the B — C delay which is 10% larger than for 
(Al + A2) — C, with a positive correlation coefficient of approximately 0.5. This implies 
that B — (Al + A2) delay is correlated with the B — C delay and anticorrelated with the 
(Al + A2) — C delay at roughly the same level. The minimum value of x 2 was 146, with 
v = 86 degrees of freedom. Scaling errors to give x 2 l v = 1 gives errors of roughly ~ 0.8 
days. 

Monte Carlo simulations of the B and C components spanning a broad range of delays 
and fitting only for the B — C delay give an rms scatter of ±3.4 days for the time delay, 
with only a weak dependence on delay length. These are considerably larger than those 
obtained from the x 2 surface, for reasons discussed by Press, Rybicki and Hewitt (1992). 
Had we simulated all three components and fit for the delays simultaneously the scatter 
might have been smaller. For the present we adopt ±3.4 days as a best estimate for all of 
the delays. 



5. Models, Distance and Hubble's Constant 

We have fit four models for the gravitational potential to the HST positions (but not 
the fluxes) of Kristian et al. (1993) for the purpose of transforming our time delay into a 
distance. All but one take the lensing galaxy to have a logarithmic 3-D potential, which 
would produce flat rotation curve for circular orbits. Model PMXS is a point mass with 
external shear. Model ISXS is an isothermal sphere with external shear. Model ISEP 
is an isothermal elliptical potential (Blandford and Kochanek 1987). The ISIS model, 
adopted by Hogg and Blandford (1995) in modelling B1422±231, uses a second isothermal 
sphere to provide shear. We identify this second isothermal with the z=0.304 group of 
galaxies (Young et al. 1981; Henry and Heasley 1986) at approximately the same redshift 
as measured for the lensing galaxy by Angonin-Willaime et al. (1993). 



The dimensionless 2-D "effective lensing potentials" ip are given by 

7 
2 



*Ppmxs(0) = & 2 lnr±^r 2 cos2(0-0 7 ) , (la) 



ipisxs(0) = br + ^r 2 cos 2(0 - 7 ) , (lb) 

fpiSEp(0) = br + 76rcos2(0 - 7 ) and (lc) 
^isis(0) = br + b'r' . (Id) 

— * 

where is angular position on the sky with respect to the lensing galaxy with polar 
coordinates r (an angle) and 0, b is the lens strength (also an angle), 7 is the dimensionless 
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shear, 7 gives the orientation of the shear (measured from W to N, consistent with 
Kochanek 1991), b' is the strength of a second isothermal sphere and r' is angular distance 
from this second deflector. The shear at any r' for model ISIS is 7 = b'/2r'. Including the 
two unknown source positions (fiw and (3n) the first three models have 5 free parameters, 
while the last model has 6 counting the P.A. and distance d of the group from the lens. 
Centering the coordinate system on the observed position of the galaxy gives 8 observable 
quantities, the coordinates of the 4 images. 

Parameter values and the source position were obtained by minimizing residuals in 
the source plane. The observed image positions were projected back into the source plane, 
as were the position uncertainties (which we took to be equal and circular). The latter 
produce source plane error ellipses with areas inversely proportional to the magnification 
which are then used as a metric in minimizing the difference between the model source 
position and the backward projections of the image positions. We take the uncertainty in 
the galaxy position to be negligible in all but one model, labelled ISIS+, for which the 
galaxy coordinates were taken as two additional free parameters. 

Results for these fits are given in Table 1. The group strengths and distances for the 
last two models are b' = 2"900 and 2'.'734 and d = U'!7 and 13"3 respectively. The strengths 
correspond to velocity dispersions for the group of 383 and 372 km/s for Q = 1. In the 
ISIS+ model the best position for the lensing galaxy was only 11 mas from the position 
reported by Kristian et al. (1993). 

The derived model parameters given in Table 1 are similar to those found by others 
(e.g. Kochanek 1991; Narasimha and Chitre 1992; Keeton, Kochanek and Seljak 1996) for 
the same system. The elliptical potential (ISEP) fit least well, and would require an E5 
or flatter mass distribution. The double isothermal (ISIS) fit significantly better than the 
isothermal with external shear (ISXS), even when the position of the "tidal" isothermal 
was fixed at the flux weighted centroid of the neighboring galaxies (c.f. Young et al. 1981) 
at P.A. -116.4° E of N and d = 1971. This corresponds to a shear at # 7 = 63.6° N of W. 
The goodness of fit, the close positional coincidence, and the plausible value of the group 
velocity dispersion are strong arguments in favor of the ISIS model. Following Narayan and 
Bartelmann (1996), time delays r are given by 



T0) = {1 + Zl) ^ 



D. 



Is 



(2) 



where (3 is the source position, the distances D are angular diameter distances, and where 
the subscripts / and s refer to lens and source, respectively. The left and right hand terms 
inside the square brackets are, respectively, the geometric and gravitational components of 
the time delay. Also shown in Table 1 are time delays in days computed for Q = 1 and 
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H = 100 km/s/Mpc, taking the galaxy redshift to be the same as the group redshift found 
by Henry and Heasley (1986), z = 0.304, and the source redshift to be 1.722. For these 
values we find 

= 30.49 days arcsec . 



c D u 

Models with (O, A) = (0.1,0.0) predict time delays longer by a factor of 1.074; those with 
(ft, A) = (0.1,0.9) predict delays longer by 1.038. 

Model PMXS predicts longer time delays because the lensing galaxy is more centrally 
concentrated. Wambsganss and Paczynski (1994) have noted this effect and emphasize 
the difficulty in constraining the degree of central concentration with only the positions of 
4 images. Model ISEP predicts longer time delays because the source must lie near the 
"diamond caustic" to produce the two brighter images. The diamond caustic for an ISEP 
model is roughly twice that of an ISXS model with the same shear, forcing the source 
further off center and increasing the difference in arrival time. A change in position of the 
lens of 11 mas produces a 3% change in the B — C time delay. Kristian et al. (1993) quote 
a 50 mas uncertainty for that position. 

Comparison of the observed B — C delay of 23.7 days with the predictions of Table 1 
gives H = 84, 44, 64, 41 and 42 km/s/Mpc for the PMXS, ISXS, ISEP, ISIS and ISIS+ 
models, respectively. The 14% uncertainty in the time delay gives a 14% uncertainty in 
each of these numbers. Interpolating between the PMXS and ISXS models, we estimate 
that taking the leading term in equations (lb-d) to vary as r 1±a2 produces an additional 
14% uncertainty in the Hubble constant. 

In all of the models the A — C delay is consistently a factor of 1.5 larger than the B — A 
delay, while the measured A — C delay is a factor of 1.5 smaller. While this discrepancy is 
only at the 2a level for our adopted uncertainties, it is nonetheless cause for concern. We 
believe that the B — C delay is less subject to systematic error than either of the shorter 
delays because of the larger number of samples per time delay. 



6. Prospects for Improvement 

There are many ways in which our time delay and distance estimate might be improved 
upon. As described above, the intranight scatter in the photometry is considerably closer 
to the photon limit than the night-to-night scatter. We expect that a new reduction 
with improved algorithms, a better profile for the lensing galaxy and perhaps with better 
flatfielding will improve the scatter and decrease the uncertainty in the time delay. Better 
algorithms may allow the inclusion of roughly 20 nights of WIYN data which fill many 
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of the gaps and have not been included in the present analysis. Finally, not all of our 
data have been reduced. Some 25 nights of data taken with the Hiltner telescope, the 
NOT, and DuPont telescope await reduction. Further monitoring would likewise reduce the 
uncertainties in the time delays and help to clarify the discrepancy between the observed 
intermediate delays and those predicted by our models. New HST data to be obtained 
during Cycle 6 is expected to reduce the uncertainty in the lens position and give an 
ellipticity for the lensing galaxy. 

We thank M. Pierce for obtaining observations of PG 11 15+080 on our behalf. PLS 
thanks G. Canalizo and E. Fulton for their early efforts in PSF fitting and lens modelling, 
and P. Steinhardt for pointing out a small numerical error in an early version of the 
manuscript. DOR thanks the Institute for Advanced Study for its hospitality. 
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Fig. 1. — An I filter direct image of PG1115+080 showing QSO components Al, A2, B, and 
C. Stars "B" and "C" of Vanderriest et al. (1986) are outside the field. 

Fig. 2. — Unshifted light curves for the PG1115+080 components components Al + A2 (*), 
C (o), and B (•) and for reference star "B" (□). 



Fig. 3. — Light curve for component C (o) shifted by 23.7 days to match that of component 
B{.). 
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Table 1: Model Parameters and Time Delays for Lens Models 
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